Electron and optical phonon temperatures in electrically biased graphene 
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We examine the intrinsic energy dissipation steps in electrically biased graphene channels. By 
combining in-situ measurements of the spontaneous optical emission with a Raman spectroscopy 
study of the graphene sample under conditions of current flow, we obtain independent information 
on the energy distribution of the electrons and phonons. The electrons and holes contributing to 
light emission are found to obey a thermal distribution, with temperatures in excess of 1500 K in the 
regime of current saturation. The zone-center optical phonons are also highly excited and are found 
to be in equilibrium with the electrons. For a given optical phonon temperature, the anharmonic 
downshift of the Raman G-mode is smaller than expected under equilibrium conditions, suggesting 
that the electrons and high-energy optical phonons are not fully equilibrated with all of the phonon 
modes. 
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The charge transport properties of the two-dimensional 
system of single-layer graphene have been the subject of 
much recent interest [I]. Graphene is characterized by 
remarkably large room-temperature carrier mobilities (> 
100 000 cm 2 V" 1 s- 1 ) i, as well as by its ability to sustain 
high current densities (> 10 s A/cm 2 ) While the 

physical processes defining the carrier mobility and scat- 
tering at low fields have been probed experimentally [7|, 
relatively little is known about the underlying physics of 
high-field transport and, in particular, the mechanisms 
responsible for the observed current saturation with in- 
creasing bias • Recent investigations indicate a crit- 
ical role for coupling of energetic carriers with the high- 
energy optical phonons [J]. The distribution of energy 
within and among the charge carriers, the strongly cou- 
pled optical phonons, and the other phonons in graphene 
has, however, not yet been established. Knowledge of 
these quantities is fundamental to establishing the ori- 
gin of current saturation. This is a subject not only of 
great inherent interest, but also of importance for emerg- 
ing applications of graphene in high-performance field- 
effect transistors [1, Hj]. Non-equilibrium behavior in the 
saturation regime, including an over-population of the 
strongly coupled optical phonon modes, has previously 
been invoked to explain the nature of high-field electrical 
transport in carbon nanotubes (Tol - [l3j . 

In this paper, we present the results of in-situ measure- 
ments of the energy distribution in key degrees of freedom 
of graphene under high electrical bias. We introduce a 
new method for analysis of the energy distribution of the 
charge carriers. The approach is based on characteriza- 
tion of light emission from the biased graphene. We find 
that over an energy range of 0.7 — 1.3 eV, the carriers 
can be described by a thermal distribution, with temper- 
atures T el reaching values as high as 2000 K before de- 
vice failure. In addition, we apply Raman spectroscopy 



to establish the population of the zone-center (G-mode) 
optical phonons T op . This temperature matches, to ex- 
perimental accuracy, that of the electrons. Based on the 
frequency shift of the G-mode phonons in the Raman 
measurements, we can also learn about the population of 
the lower-energy anharmonically coupled modes. They 
too are significantly heated, but the observed frequency 
shift is smaller than the value expected under equilibrium 
conditions, indicating that the strongly coupled optical 
phonons may not be in full equilibrium with the other 
phonon modes of the graphene lattice. This latter ob- 
servation is compatible with the trend reported in the 
recent work of Chae et al. @, 03] ■ 

In our experiments, graphene monolayers were pre- 
pared by mechanical exfoliation [l[ of kish graphite. The 
layers were deposited onto highly doped Si substrates 
covered with a 285 ± 5 nm layer of thermal Si02, used 
as a gate dielectric. Cr/Au electrical contacts were fabri- 
cated using standard electron-beam lithography followed 
by metal evaporation. Graphene channels having nearly 
rectangular shape were chosen for our devices to avoid 
any additional etching processes. The monolayer thick- 
ness of the graphene samples was confirmed by analysis 
of the Raman 2D mode [15J . The graphene devices were 
held at room temperature and under vacuum (at pres- 
sures below 10~ 7 Torr). The optical radiation emitted by 
electrically biased graphene was collected in a backscat- 
tering geometry using an achromatic objective lens (x50, 
NA=0.5). The light was then dispersed onto a charged- 
coupled device array by a single-pass optical spectrom- 
eter. The same experimental configuration was used to 
record Raman scattering spectra. In that 2.33 eV 

laser beam was focused onto a 2 /^m diameter spot at 
the center of the graphene channel. The laser power ab- 
sorbed by graphene [3j[I3| was very low (< 2 kW/cm 2 ) 
to avoid any laser-induced heating of the sample [18]. In 
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FIG. 1: (color online) (a) Source-drain current- voltage char- 
acteristics at zero gate bias of the 3.6x1.6 /xm 2 graphene 
channel shown in the top inset (the scale bar is 2 (jm). The 
lower inset shows the back-gate dependence of the channel 
conductivity (a), (b) Spectral radiance of the same graphene 
sample as a function of the electrical power per unit area P 
dissipated in the channel (symbols). The spectra are fit to 
Planck's law (solid lines). In the fits, the thickness of the 
SiC>2 layer was fixed at 280 ± 2 nm. The dips around 2.33 eV 
are artifacts arising from the notch filter used for the Raman 
measurements [Fig. [2] . (c) Radiance for photon energies from 
1.33 eV— 2.59 eV as a function of the electronic temperature 
(T e( ) deduced from the fits in (b). The solid line is a fit based 
on Planck's law. 



order to normalize the experimental data, we recorded 
the spectral response of our system using a calibrated 
tungsten-halogen lamp emitting the spectrum of a black 
body at a temperature of 3400 K. 

Figure QJ, displays the source-drain current-voltage 
(I SD — V SD ) characteristics of a typical graphene channel 
measured at zero gate bias {V G ). The minimum conduc- 
tivity point was found at V G = — 7 V, which suggests a 
residual charge density in the range 10 11 cm -2 — 10 12 cm -2 
at V G = V . In agreement with previous transport mea- 



surements |3j-|6j], we find that the current through the 
channel, I SD , shows saturation behavior in the high-bias 
regime [Fig.QJi]. During the measurements, we observed 
slight drifts in l SD , which are attributed to bias- and 
laser-induced shifts of the Fermi level. Thus, in our study 
of the electron and phonon populations, we present our 
results as a function of P, the average dissipated electri- 
cal power per unit area of the channel. 

For sufficiently high P, optical emission from graphene 
can be readily observed even in the absence of optical 
excitation. The raw optical emission spectra were nor- 
malized by the spectral response of our detection system 
and, importantly, by the computed energy-dependent en- 
hancement factor arising from multiple reflections in the 
graphene/SiC>2/Si multilayer system [16|. The corrected 
spectra are shown in figure QJb) as a function of P. We 
observe a rise of the spectrally integrated intensity by 
nearly three orders of magnitude [Fig. [Tic)], while P only 



increases by a factor of w 3. This enhancement is associ- 
ated with pronounced blue-shifts of the emission spectra 
with increasing P. We attribute the optical emission to 
bias-induced spontaneous emission from energetic charge 
carriers in the graphene channel. In order to investigate 
this phenomenon quantitatively, we consider the case of 
thermalized charge carriers and fit the corrected spectra 
using Planck's law for the spectral radiance it of a grey 
body, namely: 



u(E,T el ) = e(E) 
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In this expression, h is Planck's constant, c is the speed 
of light in vacuum, E is the energy of the emitted photon, 
and T, 



el is the temperature of the electronic system. Fol- 

we assume 



17] 



lowing recent absorption measurements 
that the emissivity e(E) of a graphene monolayer is con- 
stant over the energy range studied here (e « 2.3% fl7l]). 
We note that since optical transitions in the visible spec- 
tral range correspond to electronic excitations in the sys- 
tem, we associate the inferred temperature with that of 
the electronic system of graphene. 

The agreement between our data and this model, 
which uses no adjustable parameter, is remarkably good 
[Fig. Qlb)]. The extracted T eZ range from 1145 K 
at the lower end of our detection sensitivity up to 
2040 K for currents near the point of device failure. 
Also, the radiance integrated for photon energies from 
1.33 eV— 2.59 eV scales well with the corresponding inte- 
gration of Planck's law [Fig.[ljc)]. Optical measurements 
over this energy range reflect the behavior of electron and 
hole populations in graphene each with half the photon 
energy, i.e., from ~ 0.7 — 1.3 eV. Our results demon- 
strate that under high electrical bias these electrons and 
holes thermalize rapidly and can attain very high tem- 
peratures. As expected for the relevant energy range, 
the electron and hole distributions do not exhibit any 
observable sign of their drift velocity from the bias field. 

We now examine the coupling between the electronic 
system and high-energy optical phonons. As shown 
in Fig. [2] the Raman spectrum of electrically biased 
graphene is markedly different from a reference spec- 



trum taken at V", 



V. When P increases, the ther- 



mal emission from the electrons produces an increasing 
background in the Raman spectra, which eventually over- 
comes the Raman signal for P > 500 kw/cm 2 . We ob- 
serve that the integrated intensity of the anti-Stokes G- 
mode increases strongly relative to that of the Stokes 
counterpart. This effect is accompanied with a downshift 
and a broadening of both G- and 2D- Raman features. 
Finally, the integrated intensity of the defect-related D- 
mode is less than 5% of that of the G-mode [Fig. [2][a) , 
inset], indicating that the high structural quality of our 
samples is preserved. We did not observe any signifi- 
cant enhancement of the D-mode up to near the device 
breakdown threshold, where an enhanced and broadened 
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FIG. 2: (color online) (a) Broad band Raman spectra of the 
same graphene channel as in Fig. [T] under zero- and high- 
source-drain biases. The inset is a close-up showing the weak 
D-mode. (b) Anti-Stokes and (c) Stokes Raman G-modes 
(symbols). The solid lines are Voigt fits. The spectra are 
vertically offset for clarity. All spectra were acquired with 
the same integration time. The dissipated electrical power 
per unit area (P) and the optical phonon temperature {T op ) 
are indicated. 



D-mode ultimately develops. This trend was observed on 
several graphene channels, with similar breakdown cur- 
rent densities of ~ 4xl0 8 A/cm 2 . 

The absolute phonon population for the zone-center 
optical phonons can be derived from our measurements 
of the Raman G-mode. Indeed, the Stokes signal is pro- 
portional to 1 + n G (T op ), where n G is the G-mode pop- 
ulation with an effective temperature T op following from 
the Bosc-Einstein distribution, while the anti-Stokcs sig- 
nal scales as n G (T op ). Thus the measured integrated 
intensity ratio of the anti-Stokes and Stokes features 
[Fig. [UJb,c)] can be written as 



I. 
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C x exp 
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where fouj G is the G-phonon energy («195meV), k B is the 
Boltzmann constant and C, is a numerical factor. 

Our physical conclusions depend critically on an accu- 
rate measurement of C 14]. Thus, we measured I a g/Ig 



at room temperature on an unbiased sample, using a very 
low laser power. We obtained a value of C = 0.85 ± 0.15, 
which we have used to extract T op . The same experimen- 
tal factor of C = 0.85 is also obtained if we assume a E A 
dependence of the Raman signal on the scattered photon 
energy [l|| and properly take into account the distinct 
spectral responses of our setup and interference enhance- 
ments at the energies of the Stokes and anti-Stokcs scat- 
tered photons [iff 

The values of T op can also be derived directly from the 
intensity I s of the Stokes Raman signal. Since the opti- 
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FIG. 3: (color online) Electronic (T e ;, squares) and optical 
phonon (T op , circles) temperatures as a function of the dissi- 
pated electrical power. The dashed line is a guide to the eye, 
based on a scaling of the temperature as \/P. 



cal phonon population n <C 1 in the absence of current 
flow, we obtain I s {T op )/I° s = 1 + n G (T op ), where is 
the Stokes intensity for V sd — 0. The extracted values of 
T op using this approach are in good agreement with those 
determined from the ratio I a g/Is- I n the following, T op 
will be taken as the mean of these two independent de- 
terminations. We find values of T op from 440 K up to 
1570 K as P increases from 40 kW/cm 2 to 500 kW/cm 2 
[Fig. land Fig.|3]. 

Our results allow us to perform a direct comparison 
of T el and T op at a function of P. Within experimental 
accuracy, we find that these temperatures are identical 
[Fig. [3]. At high P, the increase in temperature with 
increasing power saturates. This behavior cannot be ex- 
plained by losses through the observed radiation. Even at 
the highest temperatures measured, the radiated power 
given by the Stefan-Boltzman law remains a negligible 
fraction (~ 4.10 -6 ) of P. The observed saturation pre- 
sumably reflects the rate of perpendicular heat flow to 
the underlying substrate [a|. 

Is this subset of high-energy optical phonons in equilib- 
rium with the lower-energy phonon modes in graphene? 
This question can be addressed by studying the spectral 
line shape of the Raman G-mode under electrical bias 
[Fig. |UJa,b)]. The Raman features are known to undergo 
a frequency downshift with increasing temperature as a 
result of anharmonic coupling [2lj |. Under equilibrium 
heating, the G-mode downshift reflects primarily the an- 
harmonic temperature of these secondary phonons into 
which G-mode phonons decay. 

For moderate temperatures, the downshift of the G- 
mode may be approximated as a linear function of the 
equilibrium temperature of the system (l8| . Such an 
approximation becomes questionable in the saturation 
regime where high temperatures may be reached Q . For 
our analysis, we have therefore made use of the theo- 
retical predictions from Bonini et al. 2l[. Although the 



theory has not been tested experimentally at high tem- 
peratures, the results are compatible with available low- 
temperature data (100 — 400 K) [l8(. We have used the 
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FIG. 4: (color online) (a) Frequency (u)q) and (b) Lorentzian 
full width at half maximum (r G ) of the G-mode feature as 
a function of T op . Data are extracted from Voigt fits, taking 
into account a constant Gaussian contribution (13 cm -1 ) due 
to our setup spectral resolution. The thick solid lines in (a,b) 
are the theoretically predicted temperature dependences of 
ojq and r G under equilibrium conditions [2lj]. The dashed 
line in (a) is a linear extrapolation of the results in ref. [2lj . 
neglecting higher-order anharmonic terms. To compute the 
solid line in (b), we considered a residual charge density of 
10 12 cm -2 and an anharmonic decay rate of 2.5 cm -1 at T = 

300 K m. 



theoretical results to estimate the downshift of lu g eq ex- 
pected if the system were in full thermal equilibrium at 
a temperature T op . As shown in Fig. 2] (a), the pre- 
dicted lj g eq are somewhat lower than the measured val- 
ues of ujq. This indicates that the strongly coupled hot 
electrons and phonons are not in full equilibrium with 
the lower energy phonons. For the highest tempera- 
ture regime(T > 900 K), we have simply extrapolated 
the results of Bonini et al.. This procedure almost cer- 
tainly underestimates the rate of downshift, since it ne- 
glects contributions from higher-order anharmonic terms, 
which would imply a larger departure from equilibrium. 
Nevertheless, we should stress that the observation of 
frequency downshifts as large as 25 cm -1 for u> G implies 
that the secondary phonon modes are significantly heated 
and may reach temperatures on the order of 1000 K. We 
may infer that heating of all accessible degrees of freedom 
in the graphene sample is significant, although complete 
thermal equilibrium does not appear to be attained. 

A further interesting point concerns the bias-induced 
broadening of the Raman features. This effect is com- 
pletely reversible and can lead to an increase by as much 
as a factor of ~ 5 for the width of the G-mode (r G ) 
[Fig. |4] (b)]. This width arises from two intrinsic compo- 
nents: coupling between G-mode phonons and resonant 
electron- hole pairs (i.e., Landau damping), and from the 
anharmonic coupling described above [21] . The first con- 
tribution would actually imply a decrease in r G for high 
temperatures. The anharmonic coupling would lead to 
broadening of the G-mode with increasing temperature. 
However, as shown in Fig. [4] (b), the magnitude of the 
broadening far exceeds that expected from the lifetime of 
the lowest-order anharmonic processes 2l|. We believe 



that the observed broadening may result from enhance- 
ment of the quasi-elastic scattering rate, as well as pos- 
sible contributions from higher-order anharmonic terms 



that become significant at high temperatures. 

Having examined the behavior of zone-center (G- 
mode) optical phonons, we finally address that of 
the zone-edge optical phonons that give rise to the 
(symmetry-forbidden) D-mode and the two-phonon 
(symmetry-allowed) 2D- mode 22, [23|. Analysis of the 



D-mode is difficult since the signal is weak and strongly 
influenced by the defect density, which may in itself be 
altered by temperature. Extracting phonon populations 
from the 2D-mode is also challenging. The 2D-mode in- 
volves electronically resonant intervalley scattering with 
large-momentum phonons [l^HjJ. Thus, as opposed to 
that of the G-mode, its intensity will reflect the com- 
petition between electron-phonon and electron-electron 
scattering processes. According to the model of [23J, we 
have I 2D tx (1 +J e l/jph) 2 (1 + n u) 2 ! where n D is the 
zone-edge optical phonon population at a temperature 
T op = T el , and j p h and j e i are, respectively, the electron- 
phonon and electron-electron scattering rates. A strong 
increase in j el is presumably the origin of the observed 
decrease of I 2D with increasing optical phonon tempera- 
ture [Fig. 13(a)]. 

In conclusion we have performed an in situ study of 
the electron and and phonon populations in graphene 
field-effect transistors operating in the high-bias regime. 
Electrons and holes with energies from 0.7 — 1.3 eV are 
found to be in a thermal equilibrium distribution and 
to achieve temperatures as high as 2000 K without de- 
vice failure. Direct measurement of the zone-center opti- 
cal phonon population shows that they are in thermal 
equilibrium with the electronic excitations of the sys- 
tem. The lower-energy phonons are also substantially 
heated, but appear to be characterized by a somewhat 
lower temperature than that of the electrons and op- 
tical phonons. This strong coupling between electrons 
and optical phonons is a hallmark of graphitic materials 
and is assuredly the main intrinsic mechanism responsi- 
ble for current saturation in graphene field-effect transis- 
tors. As regards the possible role of Si02 surface phonons 
in the microscopic processes underlying current satura- 
tion , our results demonstrate that, although 
extrinsic contributions from such modes cannot be ruled 
out, large populations of optical phonons (up ton~ 0.5) 
are present in the high-bias regime. 
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